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A review of the modern theory of hadronic atoms based on the non-relativistic 
effective Lagrangian approach and ChPT, is given. As an application of the general 
framework, we consider the decay of the 7r+7r _ atom into two neutral pions, and 
the energy-level shift of the ir~p atom. The analysis of the decay width of the 
7T+7T - atom has been carried out at the first non-leading order in isospin breaking, 
and at 0(e 2 p 2 ) in ChPT, that results in an accurate theoretical prediction for 
this quantity. The study of the ir~p atom energy-level shift within the same 
approach clearly demonstrates the necessity to critically reaccess the values of the 
ttN scattering lengths, extracted from the energy-level shift measurement by means 
of the potential model-based theoretical analysis. The construction of short-range 
hadronic potentials from field theory is also discussed. 

1 Introduction 

Recent years have seen a growing interest in the study of hadronic bound 
states - so called hadronic atoms - that provides an extremely valuable piece 
of information on the behavior of QCD at a very low energy. At CERN, 
the DIRAC collaboration El aims to measure the 7r + 7r~ atom lifetime to 10% 
accuracy. This would allow one to determine the difference ao — a2 of tttt 
scattering lengths with 5% precision. This measurement provides a crucial test 
for the large/small condensate scenario in QCD: should it turn out that the 
quantity ao — a* is different from the value predicted in standard ChPTa, one 
has to conclude! that spontaneous chiral symmetry breaking in QCD proceeds 
differently from the widely accepted picture. In the experiment performed at 
PSI cru, one has measured the strong energy-level shift and the total decay 
width of the Is state of pionic hydrogen, as well as the Is shift of pionic 
deuterium. These measurements yield isospin symmetric ttN scattering lengths 
to an accuracy which is unique for hadron physics. A new experiment on pionic 
hydrogen at PSI has recently been approved. It will allow one to measure 
the decay A^- p — ■> ir°n to much higher accuracy and thus enable one, in 
principle, to determine the ttN scattering lengths from data on pionic hydrogen 
alone. This might vastly reduce the model-dependent uncertainties that come 
from the analysis of the three-body problem in A 7[ - d . Finally, the DEAR 
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collaboration El at the DA<I>NE facility plans to measure the energy level shift 
and lifetime of the Is state in K ~p and K~d atoms - with considerably higher 
precision thandn the previous experiment carried out at KEK0 for K~ p atoms. 
It is expectedEI that this will result in a precise determination of the / = 0, 1 
S-w&ve scattering lengths. It will be a challenge for theorists to extract from 
this new information on the KN amplitude at threshold a more precise value 
of e.g. the isoscalar kaon-sigma term and of the strangeness content of the 
nucleon. 

In order to carry out the precision experimental tests of QCD mentioned 
above, on the theoretical side one faces the problem of finding the suitable 
field-theoretical framework for the description of the measured characteristics 
of hadronic atoms - energy levels and decay probabilities. In this work we shall 
report on the recent progress in this direction. The main message which will 
be delivered, is that a rigorous theory of this sort of bound states based on the 
merger of ChPT and the non-relativistic effective Lagrangian technique that 
was proposed originally by Caswell and Lepage B to study QED bound states 
in general, can be indeed constructed. This theory allows for a systematic 
expansion of the bound-state observables both in quark mass rh — ^ (m u +rrid), 
and the isospin-breaking parameters: fine structure constant a and the quark 
mass difference md — m u (for simplicity, we consider an SU(2) case here). 

Hadronic atoms appear to be loosely bound systems of hadrons that are 
formed mainly by the static Coulomb force. The Bohr radius of this sort of 
bound states is of order of a few hundreds of Fm, and the average 3-momenta 
of constituents lie in MeV range. For the above reason, it is evident that the 
non-relativistic framework provides the most natural and economical tool for 
handling such bound states. Rclativistic corrections are taken into account 
perturbatively, up to any given order in the expansion in inverse powers of 
masses. Below, we briefly outline the main ingredients of the non-relativistic 
effective Lagrangian approach to bound states, without going into details. 

• The non-relativistic Lagrangian describing interactions between hadrons 
and photons, is built from the non-relativistic hadron fields and the photon 
field. This Lagrangian consists of an infinite tower of all possible operators, 
with an increasing mass dimension - all operators allowed by discrete symme- 
tries and the gauge invariance should be included, with an a priori unknown 
couplings. In actual calculations, only a few low-dimensional ones matter - 
higher-dimensional operators contribute to higher powers in a in the bound- 
state observables. This feature of the effective theory goes under the name of 
"power counting" . Further, the Lagrangian does not include, by definition, the 
operators that change the number of heavy particles (hadrons). 

• Loop corrections to the scattering amplitudes in the non-relativistic the- 
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ory are calculated in a standard manner, by using the Feynman diagrammatic 
technique. There is, however, one important modification. It is well known 
that in the non-relativistic theory in the presence of light particles (photons) 
the Feynman integrals should be properly regularized in order to avoid the 
contribution from the loop momenta at the hard scale - otherwise, loop correc- 
tions to the Green functions would lead the the breakdown of counting rules 
in the non-relativistic theory. A suitable regularization procedure built on the 
top of the Feynman rules in the non-relativistic theory is provided by so-called 
"threshold expansion" that enables one to disentangle the contributions com- 
ing from different regions of loop momenta, by expanding the integrands - in 
the dimensional regularization - in all possible small kinematical variables. 
Next, one has to systematically remove hard- momentum contribution from 
the integrals, which at low energies is given by a polynomial in external mo- 
menta, and can be absorbed into the renormalization of the couplings in the 
non-relativistic Lagrangian. For a more detailed discussion of the problem, 
seeB. 

• Couplings in the non-relativistic effective Lagrangian are determined 
from the matching to the relativistic theory. These couplings are not necessar- 
ily real, since the physical decay processes where the number of the hadrons 
is not conserved, contribute to the imaginary parts of these couplings. The 
matching condition determines the couplings in terms of threshold parameters 
of the physical hadronic scattering amplitudes. It is crucial to stress that the 
matching condition does not imply the chiral expansion of the amplitudes - 
formally, it is valid in all orders in the chiral expansion. 

• After setting the parameters in the non-relativistic Lagrangian, we turn 
to the bound states in the theory. The Feshbach formalismlHl that allows one 
to separate the bound-state pole in the scattering amplitudes, turns out to be 
very convenient for this purpose. The real and imaginary parts of the pole 
position on the second Riemann sheet of the complex energy plane coincide, 
by definition, with the energy and the decay width of the (metastable) bound 
state. The perturbative framework for determining the pole position coincides 
with the conventional Rayleigh-Schrodinger perturbation theory. At the end, 
the characteristics of the bound state are determined through the couplings in 
the non-relativistic amplitudes. With the use of the matching condition, these 
characteristics can be further expressed via the scattering amplitudes in the 
relativistic theory. 

• Given the chiral expansion for the scattering amplitude, it is possible 
to obtain the chiral expansion for the bound-state observables which can be 
reorganized in the expansion in m and isospin-breaking parameters. The isopin 
symmetry world is, by convention, defined as the one where the masses of pions 
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an nucleons coincide with the charged ones, and the values of all low-energy 
constants remain the same. In this vein, one may achieve an unambiguous 
separation of the isospin-breaking effects and extract the isospin-symmetric 
hadronic scattering lengths directly from the hadronic atom measurements. 

After having described the general framework, we consider to the particular 
systems. 

2 Decay of the 7r + 7r atom into 7r°7r° 

Recently, using the non-relativistic effective Lagrangian framework, a general 
expression for the decay width r^ 2jr ^ 7r o 7r o of the Is state of the n + ir~ atom was 
obtained at next-to-leading order in isospin-breakinglij. Numerical analysis of 
this quantity was carried out at order 0(e 2 p 2 ) in ChPTEj. TTipsp imEstiga- 
tions have confirmed and generalized the results of earlier studies EJ't£l. The 
expression for the decay width at the first non-leading order in isospin breaking 
has the form 

2 

r^^^o = - QfVArUl + k^) , (1) 

K ™ = §m 2 ( ~ a ° + 2ct2 '' 2 " IT ^ na - 1 ^ 2a ° + + ( md _ m «) 2 ) • 

Here p* = (M 2 + - M 2 - ^M^a 2 ) 1 / 2 , and a/, (I = 0,2) denote the (dimcn- 
sionless) strong 5*-wave tttt scattering lengthaJn the channel with total isospin 
/. The quantity A^n is obtained as follows E3. One calculates the relativistic 
amplitude for the process ir + ir~ — ► tt q tt at 0(a, {rnd — m u ) 2 ) in the normal- 
ization chosen so that at O(l) the amplitude at threshold coincides with the 
difference ao — 0,2- Due to the presence of virtual photons, the amplitude is 
multiplied by an overall Coulomb phase 9 C that is removed. The real part of the 
remainder contains terms that diverge like |p| _1 and In2|p|/M 7r + at |p| — > 
(p denotes the relative 3-momentum of charged pion pairs) . The quantity Ann 
is obtained by subtracting these divergent pieces, and by then evaluating the 
remainder at p = 0. 

Re (e~^ t„„) -► p- + b 2 In |M. + J^L- + ■ ■ ■ (2) 
IPI M n + oM 7r+ 

As it is seen explicitly from Eq. (Q), one can directly extract the value of 
A n n from the measurement of the decay width, because the correction K vv is 
very small and the error introduced by it is negligible. In order to extract strong 
scattering lengths from data, one may invoke ChPT and to relate the quantities 
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A-ktx and do — a 2 order by order in the chiral expansion c3. This requires the 
evaluation of isospin-breaking corrections to the scattering amplitude. At order 
0(e 2 p 2 ) in chiral expansion we obtain 

A^. = a - a 2 + , e nn = (0.58 ±0.16) ■ 1CT 2 , K n7r = 1.15 • 10~ 2 . (3) 

The lifetime of the ir + ir~ atom is predicted to be 

r 27r o = (2.91 ±0.09) • 10" 15 s, (4) 

and the correction to the leading-order in isospin-breaking r 2;r o = (1 + 5p) 
where T^ = ^a 3 p*(a - a 2 ) 2 , equals to S r = 0.056 ± 0.012. Note that 
above we have used the values of irir scattering lengths ao = 0.220, a 2 = 
—0.0444, A(ao — a 2 ) = 0.004, obtainecLon the basis of two loop calculations in 
ChPT and dispersion relations analysisE3. For the various low-energy constants 
entering the expression for , the same values as in Ref. c3 were used. These 
values of the low-energy constants will not necessarily coincide with the ones 
that will emerge from the analysis of scattering lengths. In addition, the errors 
in the scattering lengths and the correction term at this stage are treated as 
uncorrelated. Only once the new values for the low-energy constants U from 
the dispersion analysis are available including error bars, one may refine the 
above preliminary analysis for the decay width. We also emphasize, that the 
corrections to the decay width at 0(e 2 z2^ are already very small, justifying 
the negligence of the higher-order terms 113. 



3 Energy level shift of the ir p atom 

The treatment of the ir~p atom problem0 proceeds along the lines very similar 
to those for tt + it~ case. Our investigations are aimed at the derivation of the 
general expression for the n~p atom energy-level shift in the Is state. The 
total shift is given by a sum of the electiomagnetic and strong pieces. Our 
calculations for the electromagnetic shift EH within a high accuracy yield the 
same result as given in Ref. □. The final result for the strong shift e\ s in the 
first non-leading order in isospin breaking is given in a form similar to Eq. ([!]) 

e ls = -2a 3 [i 2 c AnN (1 + K^n) (5) 

KirN = -2a^ c (lna - l)ArJV + o(a,m d ~ m u ) , 

where [l c denotes the reduced mass of the n~p pair, and the quantity A^n is 
defined analogously to A-n^. To calculate this quantity, one has to evaluate 
the n~p — > ir~p relativistic scattering amplitude at 0(a,rrid — m u ), drop all 
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diagrams that are made disconnected by cutting one photon line, and discard 
the spin-flip piece. The remainder is denoted by t^N- The regular part of t^N 
at threshold defines the quantity A^n in analogy to Eq. (^) 

Rc (e~ 2 ^ t vN ) -» + B 2 In M - — A* N + ■ ■ ■ . (6) 

IpI Mc 

In order to extract the value of the S-wave ttN scattering lengths , ajT, 
from the ir~p energy shift measurement, one may again resort to ChPT, to 
calculate the isospin-breaking corrections to the ttN scattering amplitude at 
threshold. The normalization of the relativistic amplitude is chosen so that 
at threshold At,n = a o+ + a o+ + 0(a, m<j — m u ). We have carried out these 
calculations at chiral order 0(p 2 ), where only the tree diagrams contribute, 
the result looks as follows 

A + , - , _ m p (8 Cl AM 2 - 4e 2 /! - e 2 / 2 ) 
Awv-a 0+ + «o + + ^v, ^n- 87r(mp + Mjr+)F 2 ' ( ? ) 

where Ci (/j) are the strong-/ electromagnetic) low-energy constants from the 
0(p 2 ) Lagrangian of ChPTlia. In order to perform the numerical analysis, one 
has to specify the values of these low-energy constants. The "strong" constant 
ci can be determined fijarn the fit of the elastic irN scattering amplitude at 
threshold to KA86 dataN: a = -0.925 GeV -1 . The value of the constant 
can be extracted from the proton-neutron electromagnetic mass difference E3: 
e 2 / 2 = (—0.76 ± 0.3) MeV. The determination of the constant fx from data 
is however, problematic. For this reason, in our analysis we have used order- 
of- magnitude estimate for this constant: — 1/ 2 < /1 < 1/2!- With these values 
of the low-energy constants, we obtain the isospin-breaking correction to the 
leading-order energy-level shift defined as e\ s = ef s °(l + Se), where ef^ = 
— 2a 3 /x 2 (aQ + + % + ), to be S t = (—4.8 ± 2.0) • 10 -2 . The large uncertainty is 
caused mainly by the poor knowledge of the parameter j\. In demonstration 
of the above discussion, in Fig. 1 we confront the results of the above analysis 
with those of the potential modelQ: S £ = (—2.1 ± 0.5) • 10~ 2 , using the same 
experimental input. As it is readily seen from Fig. 1, the systematic error in 
the potential approach is grossly underestimated, that already indicates at the 
necessity to critically reaccess the values of the S'-wave ttN scattering lengths 
quoted in Ref. a. In addition, it remains to be seen, how our results will be 
altered by the loop corrections in ChPT. A precise determination of the low- 
energy constant fx, using either sum rules or invoking various models, is also 
desirable. 

As it was mentioned in the Introduction, a new experiment at PSI is ap- 
proved, that is aimed at a precise determination of the ir~p atom decay width. 
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Figure 1: Determination of ttN scattering lengths from the pionic hydrogen and pionic 
deuterium measurements. Solid line corresponds to the energy-level shift calculations at 
0(p 2 ) in ChPT, and the dashed lines - to the potential model results 

At the first non-leading order in isospin breaking, the width is determined 
by the charge-exchange amplitude ir~p — > 7r°n where the low-energy constant 
fi does not appear. Consequently, the calculations for this quantity can be 
carried out with a less theoretical uncertainty. 

4 Construction of the potentials from field theory 

The study of the properties of the hadronic atoms within the potential ap- 
proach has a decades-long history. As a general rule, the predictions made 
within the potential model badly deviate from the ones obtained within the 
field-theoretical approach based on ChPT. For example, the potential model 
predicts the isospin-breaking correction 5r to the it + tt~ aiom decay width with 
an opposite sign and with the same order of magnitudeEH. The predictions for 
the ir~p atom energy-level shifta have been already considered in the previous 
section. The reason for this discrepancy is now well understood. Namely, the 
potential model does not take into account a full content of isospin-symmetry 
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breaking in QCD. In particular, the effect of the direct quark-photon inter- 
actions (encoded in the "electromagnetic" low-energy constants of ChPT), as 
well as the effect stemming from the dependence of scattering amplitudes on 
the quark masses, are not included. As these effects contribute a bulk of the 
total isospin-breaking correction term in the field-theoretical approach, it does 
not come to our surprise that the predictions of both approaches substantially 
differ. 

Since the isospin-breaking effects discussed above lead, in the language 
of the potential scattering theory, to the modification of the short-range part 
of the hadronic potential, it is evident that in order to bring the potential 
model in conformity with the field-theoretical treatment, one has to assume 
that short-range hadronic potentials contain an isospin-breaking piece. It is 
natural to seek a derivation of the potentials that are used in the potential 
model, on the basis of ChPT. In a slightly more restricted context, one may 
ask, how the isospin-breaking part of the short-range "strong" potential is 
obtained from ChPT, when the isospin-symmetric part is already known to fit 
well ChPT predictions (we recall that the isosaiiubreaking part is assumed to 
vanish identically in existing potential modelsu'Eil) . 

It is widely presumed that the potential constructed from the field theory 
will be necessarily singular in the position space and will require some kind of 
regularizationEi We argue that this is not necessarily the case: almost ^ny 
well-behaved short-range potential, including those that were used in Refs.QO, 
can be generalized to include properly the full content of isospin-breaking ef- 
fects in ChPT. 

The key observation that leads to the above conclusion, can be summarized 
in the so-called universality conjecture. This conjecture - completely in spirit 
of the low-energy effective Lagrangian approach to bound systems - states that 
the bound-state energies in the field theory, and in the potential model are the 
same at the first order in isospin breaking, provided the threshold amplitudes 
calculated in these two theories, coincide. We shall ensure the universality for 
the case of a simple one-channel model where the interaction Hamiltonian is 
given by a sum of Coulomb and short-range interactions Hi = — qt _1 + U. 
The short-range potential U, in general, contains isospin-breaking effects. 

With a given interaction potential, one may evaluate the energy-level shift 
of the ground state of the bound system. The equation for the position of the 
bound-state pole in the (complex) energy plane is given bytHI 



z-E -(* \t(z)\Vo)=0, 

t(z) = U + U(z - Ho)- 1 ^ - l*o>(*o|)r(z) , 



(8) 
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where Eq denotes the unperturbed Coulomb energy of the ground state, and 
\I>o stands for the unperturbed wave function. The iterative solution of Eq. (||) 
yields usual Rayleigh-Schrodinger perturbation series. Further, it turns out 
useful to introduce the scattering matrix on the short-range potential only 



t(*)=U + U : — t(z), to EE <q|t(z)|p) 

Z — tin 



(9) 

p=q=0, z=0 



Here to denotes the value of the scattering amplitude at threshold (for conve- 
nience, we have shifted the threshold to z = 0). 

The perturbative solution of the equation for the bound-state energy at 
0(a 4 ) gives 

en = ^ (to - ^ tg (l + In -L.) + fi [t] + ^ m <>) } ■ (10) 

where Q and 1Z are certain known functionals of the scattering matrix t. We 
do not display the explicit expressions here. The dependence on the arbitrary 
cutoff parameter b in the functional 1Z cancels with the logarithmic dependence 
in the second term, so that de\ s /db = 0. 

At the next step, we consider the full scattering amplitude defined through 
the Lippmann-Schwinger equation T(z) = Hi+Hi(z— Hq) _1 T(2). The thresh- 
old behavior of T(z) on energy shell is given by (cf with Eq. ra)) 



Re[e-^<p|(T(*)-H,)|q>] 



|p|=|q|i z=p 2 //Jc 



5i + s 2 inM 

P Mc 



2-7T 

-—A + 0{\p\) + o(a) (11) 
A = to - t g ln J_ + rr^ Q[ t ] + ^[ t; 6] + (a) , (12) 



/jL c IT i[l c IT 7T J 



with the same functionals Q and 1Z. 

If one now expresses the energy level shift in terms of the threshold ampli- 
tude A, one arrives at exactly the same expression (||) as in the field-theoretical 
framework - in accordance with the universality conjecture 



e ls = -2a 6 t/ c A (1 - 2a/A c (ln a - 1)A) + o(a 4 ) . (13) 

Based on the universality conjecture, we can provide a constructive algo- 
rithm for the derivation of the isospin-breaking part of the short-range po- 
tential U from ChPT. The amplitude at threshold in the latter is generally 
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given by A = Ao + A\ + ■ ■ •, where Aq, A\ denote the isospin-conserving 
(breaking) parts of the amplitude, and ellipses stand for higher-order terms 
in isospin breaking. In order to ensure the inclusion of the full content of 
isospin-symmetry breaking in ChPT into the potential model, it thus suffices 
to match the amplitude A in both theories. The problem evidently has too 
much degrees of freedom. The short-range potential is also given by the sum 
of isospin-conserving and isospin-breaking pieces U = Uo + Ui + • • •, and for 
our purposes the following ansatz is sufficient: U = (1 + A)Uo, where the sole 
coupling A will be determined from matching the isospin-breaking pieces of the 
threshold amplitude. 

If one defines the scattering amplitude in the limit of no isospin breaking 
through t(z) — Uo+Uo(z— Ho) _1 t(z), The matching condition may be written 
recursively, in terms of t 

Aq = t , (14) 

-*±A X = A(t + S[i]) - tV In ± + ^ Q[i] + ^ K[i: b] . 
where S again stands for a certain known functional. The matching condi- 



tion (14) solves our problem completely - the bound-state energies calculated 
with the use of the "corrected" potential coincide, by definition, with those 
calculated on the basis of ChPT. 



5 Conclusions 

The approach based on non-relativistic effective Lagrangian technique and 
ChPT, provides one a powerful tool to systematically calculate the charac- 
teristics of loosely bound states of hadrons. With the use of this approach, the 
7r + 7r~ atom decay problem is now completely understood, both conceptually 
and numerically. Certain theoretical effort will be still needed to extract the 
precise values of ttN scattering lengths from past and future measurements at 
PSI. The treatment of other systems, like kaonic atoms that will be measured 
by the DEAR experiment at DA$NE, is foreseen within the same framework. 

Further, within the present approach, one may establish a constructive al- 
gorithm for the derivation of the short-range hadronic potentials from ChPT. 
The algorithm is based on the universality conjecture that has been discussed 
above, for a simple one-channel case. We hope that - after the suitable gener- 
alization - the approach based on the universality might be also useful for the 
analysis of irN scattering data near threshold, in what concerns the study of 
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the isospin-breaking effects in the ttN amplitude. 
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